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Abstract 



For the purpose of understanding the quantum behavior such as quantum decoherence, fluctuations, 



dissipation, entanglement and teleportation of a mesoscopic or macroscopic object interacting with a general 
environment, we derive here a set of exact master equations for the reduced density matrix of N interacting 

en- 
harmonic oscillators in a heat bath with arbitrary spectral density and temperature. Two classes of problems 



of interest to us which these equations can be usefully applied to are that of the quantum dynamics of 
nanoelectromechanical oscillators and the entanglement evolution of multi-partite macroscopic states such 
as quantum superposition of mirrors in a high Q cavity. To address a key conceptual issue for macroscopic 
quantum phenomena we examine the conditions for an assumption often implicitly made in these studies to 
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be valid, namely, that the quantum behavior of a macroscopic object in an environment can be accurately 
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can be used to calculate the uncertainty principle governing a macroscopic object at finite temperature. 

O 



represented by only treating the dynamics of its center-of-mass variable. We also mention how these results 
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I. INTRODUCTION 

In an earlier paper [l| we showed the derivation of an exact master equation for two coupled 
quantum harmonic oscillators interacting bilinearly with a common environment made up of n 
harmonic oscillators at an arbitrary temperature for a general spectral density function. (This is 
referred to as a 'general environment' in [2].) This equation can be applied to the analysis of model 
problems in macroscopic quantum phenomena (MQP) involving two harmonic oscillators, either 
mechanical such as the superposition of two mirrors 3|, y, |5], or one of them of electromagnetic or 
superconducting flux origin, as in nanoelectromechanical (NEM) resonators [6] or SQUID-resonator 

In this paper we take a small step towards treating the quantum properties of a mesoscopic 
or macroscopic object in a general environment by providing the technical base for such studies. 
We consider a system modeled by N harmonic oscillators (NHO) interacting with a heat bath 
consisting of n harmonic oscillators (HOB). The aim is to delineate the conditions upon which 
the mechanical and statistical mechanical properties of this quantum object can be described in 
terms of its center-of-mass (COM) variable by a master equation for the reduced density matrix 
with the bath variables integrated out, and for such conditions, derive an exact master equation 
for a bath with arbitrary spectral density and temperature. The motivation for this work has both 
conceptual and practical underpinnings. At the conceptual level we want to examine the validity of 
an implicit assumption made in many MQP investigations, namely, that the quantum mechanical 
behavior of a macroscopic object like the NEM or a C60 molecule [§[ Q], placed in interaction 
with an environment, behavior such as quantum decoherence 10], fluctuations and dissipation, 



entanglement and teleportation, can be captured by its COM behavior. For convenience we refer 
to this as the 'COM axiom'. This assertion is intuitively reasonable, as one might expect it to be 
true from normal- mode decompositions familiar in classical mechanics, but when particles (NHO) 
interact with each other (such as in a quantum bound state problem) in addition to interacting 
with their common environment, all expressed in terms of the reduced density matrix, it is not 
such a clear-cut result. At least we have not seen a proof of it 3C|. 



Another important issue is the demarkation between microscopic, mesoscopic and macroscopic. 
Assuming that the object is made up of N physically relevant quantum particles (e.g. atoms, 
forgetting about the tighter-bound substructures), starting with N = 1 which we refer to as micro- 
scopic, the question is: At what number of N will one begin to describe the object as mesoscopic 
with a qualitatively distinct behavior from microscopic, and likewise for macroscopic? In clas- 
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sical statistical mechanics, this issue underlies the important attempt to derive from molecular 
Hamiltonian dynamics (with deterministic chaos) the thermodynamic and kinetic properties of a 
gas (with N molecules), such as transport functions, and their dynamics, which possesses salient 
dissipative and time-asymmetric features. To provide a quantitative analysis of such issues one 
needs to work with a stochastic equation for the N particles so as to be able to see the cross-over 
behavior between any two of these three regimes as one varies N, and its thermodynamic behavior 
as one varies the temperature T. In this paper we derive such an equation for NHO in HOB, but 
will leave the analysis of this theoretical issue to later investigations 12 ] . 



On the applied side, the master equation for NHO in a HOB is useful for a range of problems 
which are of experimental interests. Detection of small displacements of a NEM resonator by 
a superconducting single electron transistor (SSET) or a biased quantum point contact (BQPC) 
is a useful scheme for probing the appearance of quantum properties of oscillators, such as 



decoherence and entanglement 



14| . noise and fluctuations [lj]], the standard quantum 



and the uncertainty principle for quantum open systems at finite temperatures 



been shown that in certain well-accessible regimes both the SSET and BQPC 



17, 
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m. It has 
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devices behave like a thermal bath (albeit in some regimes it shows strong back-action effects which 



have been suggested as a creative way to cool the resonator 
equations have been derived for these systems in 
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22J). Master and Fokker-Planck 



23]. 



Another interesting setup is a linear array of NEM resonators. Studying the entanglement 
transport in chains of mechanical oscillators, Eisert et. al. 24| showed that the motions of distant 



oscillators can be entangled without the need for control of individual oscillators and without 
assuming any direct interaction between them. Our master equation for N harmonic oscillators in 
a general environment can be usefully applied to this problem to check on the distinct robustness 
of entanglement in the canonical coordinates found by these authors. 

Technically, our derivation of the master equation for the NHO system is similar to the simple 



method we used in 



for the case of two coupled harmonic oscillators in a HOB. The task rests on 
finding a suitable canonical transformation which preserves the structure of the Poisson brackets. 
In Sec. 2 we define our model and outline the procedure. In Sec. 3 we write down the evolution of 
the density matrix and derive the master equation which has the same form as the Hu-Paz-Zhang 



(HPZ) equation 



2a. 1261. |27J. In Sec. 4, we derive the corresponding Fokker-Planck equation. In 
Sec. 5 we discuss the general features of these equations. The overall characteristics of the results 
derived in this paper may posteriori be deduced from intuitive reasoning. However, the fine points 
present in the full solution are not so easily obtained. They are needed to address the validity of 
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the COM axiom we posed above, i.e., under what conditions can one presume that the center of 
mass coordinate of a macroscopic object is the one most sensitive to the environmental influence? 
We discuss this issue in the last section. In the Appendices we derive the canonical transformation 
for a general NHO system, and show the explicit construction procedures for the cases of N = 2 — 5. 



II. THE MODEL 



For the investigation of environmental influences on the quantum properties of a macroscopic 
object we consider the quantum Browninan motion (QBM) of N identical harmonic oscillators 



2m- 



with mutual interactions coupled to a collection of n oscillators making up the environment 
The generalized QBM Hamiltonian is: 

TV N 
H sys = ^2(-Mx* + -Mtfxt)+ \),Lr, ;■/) (1) 

n 1 1 

#bath = ^2 m ^ + 2™^ 2) (2) 
i=l 

N n 

Him = ^ E Ci 3 ( Xiq y> ( 3 ) 
i=l j=l 

For simplicity, let us assume that = Cj,\/i = 1,...,N i.e. all the system harmonic oscillators 
couple to the bath with equal strength. Hence the interacting Hamiltonian can be rewritten more 
compactly as 

N n N n 

i=l j=l i=l j=l 

With this generalized N harmonic oscillator (NHO) QBM model, we can show the following: 

(1) For any given finite iV sets of canonical coordinates (xi,Pi),i = 1,...,N, we will give a 
procedure (or algorithm) to construct another set of coordinates (Xj,Pj),i = 1,...,N where X\ = 
jf(x\ + ... + xjy),Pi = P\ + ... + Pn are the center-of-mass coordinate and total momentum of 
the system respectively. This transformation matrix T which transforms x to X, i.e. Xj = TijXj, 
satisfies |detT| = 1. The choice of {Xi} is not unique. However, our method gives a systematic 
way of construction and is thus quite useful. We also provide the transformation for the effective 
masses. The explicit construction is shown in Appendix A. 
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(2) The quadratic part of the system Hamiltonian can be shown to transform as 




where Mi = Yli=i M = NM is the total mass of the system. 

Under this transformation, the Poisson brackets and hence the commutation relations are pre- 
served, 

[xi, Pj] = [Xi, Pj] = ih5ij, [ Xi , xj] = [Pi, Pj] = [Xi,Xj] = [P^Pj] = 0. (6) 

Hence our transformation is a canonical one. Most importantly, we will show that the interaction 
potential between the N harmonic oscillators in the new coordinate system is independent of X±, 
the center-of-mass coordinate. 



Lemma: If the potentials Vij among the system oscillators Xi and Xj are functions of X{ — Xj 
only, then = 0. Hence the total potential Y^Hvi+i Viji x i ~ x j) is independent of X±. The proof 
of the lemma will be given in Appendix B. 

Combining all the above properties, we can rewrite the original Hamiltonian as follows: 

N 1 1 N 

H sys = Y.^2 M ^ + 2 M ^ } + E V n( x i- X j) ( 7 ) 
i=i " ij=i t i^j 

N P 2 1 

= ^{^ + 1 M i n 2 xf) + v(x 2 ,...,x N ) (8) 

1=1 l 
62 i N 62 i 

= ^ + -M^Xf + Y J (^ + ^ 2 Xf) + V(X 2 ,...,X N ) (9) 
= Hi + H 2 , (10) 



where 



H^^Ly-M^Xl (11) 



N P 2 1 

H2 = E(d=r + o M ^ 2x ?) + V ^ •••> ( 12 ) 

i=2 2Mi 



N n N n n 

i=l j=l i=l j'=l j'=l 
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Note that the heat bath has the spectral density J (to): 



n £j 2 



(14) 



j=l 3 3 

which differs from the original heat bath density by a numerical factor N 2 . 
The total Hamiltonian can then be written as 

-fftot = -f^sys + #bath + -f^int 

= H\ + H2 + -ffbath + #int 
= H cm + H' + #bath- 

H cm = H 1 + H mt , H' = H 2 . Note that [H', H cm ] = [H', H hath ] = 0. 

III. DENSITY MATRIX AND MASTER EQUATION 

In our derivation, we shall make the following two assumptions: (1) The system and the envi- 
ronment are initially uncorrelated. (2) The heat bath is initially in a thermal equilibrium state at 
temperature T = {ks(i)~ l ■ 



A. The density matrix 

The density matrix for the total system develops in time under the unitary evolutionary oper- 
ator: 



p(t) = exp 
= exp 
= exp 



.Htatt_ 



p(0) exp 



.Htott 



h 



. {H cm + H' + -ffbath)* 



h 



p(0) exp 



. (H cm + H' + #bath)£ 



H't 



exp 



. (^cm + #bath)^ 



P(0) exp 



n 



■ {H CIn + ffba th)^ 

; h 





r H't] 


exp 





(15) 



In the third equality, we have to use the condition that [H cm , H'] = [i?bath, H'] = 0. (Note 

that for general operators A, B, we have the Baker-Campbell-Hausdorff formula: e A e B = 
e A+B+\[A,B\+^{[A\A,B\]+[B,[B,A\])+...^ 

If we define 



p(t) = exp 



— 1- 



■ {Hc m + H ha-th )t 

h 



p(0) exp 



. (^cm + ^bath)^ 

h 



(16) 
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then 



pit) = exp 



.H't 



p(t) exp 



.H't 



(17) 



p r = Tr hat hP(t) = exp 



.H't 



Tr hath p(t) exp 



(18) 



In a similar manner as in [1], we can get the exact master equation for the arbitrary N oscillators. 



B. The Master equation 



Tracing over the heat bath leads us to a HPZ type master equation [2| for Xi,P±: 

Pr = ^[Hl,Pr] + ^§-[Xl~Pr] + ^r[X U {P 1 ,p r }] + ^-[X U [P^Pr]] ~ ^[^l, [X X ,p r ]\ (19) 



ih 



2ih 



2ih 1 



h 2 



where H\ is the Hamiltonian for Xi,P± only. This has the same form as the HPZ equation for 
the QBM of a single harmonic oscillator X\,P\ interacting with a general heat bath. Note that 



coefficients satisfy the same equations as listed in Refs 2j or 
constants and masses. 



2q |. but with different coupling 



From the evolution equation (fT7|) . the required master equation for p r (t) is obtained, 

Pr = ^[H sys , Pr ] + ^[XlPr] + {Pl,Pr}] + [h, Pr]] ~ Pr]] (20) 

The only difference between Eq. (|2U|) and Eq. (|19|) is that the unitary evolution is modified by the 
remaining N — 1 fictitious harmonic oscillators (Xj, Pj),j = 2, 3, N. 
In terms of the original variables x%, ...,xjsr, Pi, Pn, we get 

Pr = T^t^y^^] + 2AT2^ [( Xl + ■•• + x n) 2 ,Pt] + ^pfc^i + - + XN,{ p i + ■■• + Pn,Pt}] 

+ n¥^ Xi + "" + XNl + "' + P/v ' pr ^ ~ J^h^ Xl + "' + XN ' ^ Xl + "' + XN ' pr ^ ^ 

This exact master equation for the N coupled harmonic oscillators in a general environment is the 
main result of this paper. 

IV. FOKKER-PLANCK EQUATION 

In this section, we present an alternative, but useful form of the master equation derived in the 
last section. We also provide the explicit expressions for the coefficients appearing in (j20j) . 



A. Fokker-Planck equation 



In terms of the Wigner function, the above master equation takes the form: 
dW J^f Pi dW . ~ dw\ 



dt ^ M dxi dPi ) 

+ un mxi + ... + w) + ... + * + 2rW ( J_ + ... + _JL) [(Pl + ... + P K )^] 
9 <9 \ 2 ~ . , , f d d \ f d d 



Note that the Wigner function is related to the reduced density matrix in the following way: 

W{ Xl ,..,x N ,P u ..,P N ,t) (23) 
= WWI dUl - dUN ^ lPl+ -" WV)/ V (*i - %,.,x N - f ; X1 + ^,..,x N + f,t) 



where we identify in Eq. f)22[) : 



a(t) = MAO(t), (24) 

&(*) = 2r(t), (25) 

c(t) = A(t), (26) 

d(t) = E(t) (27) 

By using this transformation, the master equation can be easily obtained: 

< d » - * * 92 92 ) Pr+ 1 -M^> 

i a 



+ ^M5fl 2 (t)(xi + ... + x N - yi - ... - y N )^(xi + ... + x N + y 1 + ... + 2/jv)Pr 

1 d d d d 

-ihT(t)(xi + ... + XN ~ VI - - - VN)7;(a 1" ••• + a a ••• - a )Pr 

2 ox i oxn oy\ oyN 

-iMT,(t)(xi + ... + x N - yi - ... - yu fpr 

+hA(t)( Xl + ... + x N - y x - ... - y N ){^- + ... + -7^— + + ... + ^-)Pr (28) 

dxi oxjy dyi dy N 

B. Coefficients 

The coefficients a(t), c(t), d(t) appearing in (|2ip or (|22p can be constructed in terms of the 
elementary functions Ui(s),i = 1, 2. The Ui(s) are defined as the functions that satisfy the following 
homogeneous integro-differential equation 

M 2 f' s 

S(s) + ft 2 £(s) + — / dAr/(s - A)E(A) = (29) 

M JO 



with the boundary conditions: 

Ul ( s = 0) = 1 , m(s = t) = 0, (30) 

and 

m 2 ( s = 0) = 0, U2 (s = t) = l. (31) 

where 

rj(s) = — cL)I(uj) sin(u;s) (32) 
Jo 

is the dissipation kernel and I{ui) = -^J(u)) is the spectral density of the environment. Note that 



the numerical pre-factor before the integration in this equation is different from that defined in 



HPZ 



a. 



This is the main difference induced by the presence of ./V harmonic oscillators. 
Let G\(s,t) be the Green function which satisfies the following equation: 
d 2 N 2 f' s 

—G 1 (s,T) + n 2 G 1 (s,T) + — J dr V (s-T)G 1 (s,T) = 5(s-r), (33) 
where G\(s,t) as a function of s satisfies the following initial conditions: 

G 1 (a = > r)=0, ^(s, r)\ s=0 = . (34) 
The Green function G2(s,r) is defined analogously. The coefficients can then be written as 

a(t) = N /' d,r,(t - s) L 2 (s) - M ^' ) ') , (35) 

«*) = £ (38) 
c(t) = [ d\Gi(t, \)u(t - A) 



o 



NM 
N 2 H r> 



M 2 



[ ds [ dr f dXr ) (t-s)G 1 (t,X)G 2 (s,T)u(T-X), (37) 

Jo Js Jo 



and 



d(t) 



M 

where v{s) is defined as 



^ j\xG> x {t,X>{t-X) 

I ds I dr [ dXr 1 (t-s)G' 1 (t,X)G 2 (s,T)u(T-X). (38) 

Jo Js Jo 



r+oo ^ 

u(s) = / du)I(u)) coth.(—hjj (3) cos(o;s) . (39) 
Jo 2 

which is the noise kernel of the environment. Here a prime denotes taking the derivative with 
respect to the first variable of Gi(s,t). 



V. DISCUSSIONS 



We end with a few technical remarks followed by two conceptual points, one referring to the 
COM axiom for the quantum dynamics of macroscopic objects, and the other, to the generalized 
standard quantum limit. 

A. Technical Remarks 

First, note that although we can show [Xi,H2] = [A, #2] = and hence [H\,H2] = 0, the 
(N — 1) set of variables Xj,j = 2, ...,7V will in general not commute with H2 because of the 
potential V(X2, -Xjv)- 

Here we assume that all the N particles are of the same mass M and have the same eigen- 
frequency Q. They both couple to the environment with equal strength. So these N harmonic 
oscillators are "identical" particles. If the N system oscillators have different masses, this becomes 
more involved. 

In the proof that the potential is independent of X\, we made no assumption about the func- 
tional form of the potential. All that was needed was that the potentials are functions of Xj — Xj 
only. This is a reasonable assumption valid for many interesting physical situations. Hence the 
range of applicability of our result is by no means overly restrictive. 

B. The COM Axiom for quantum dynamics of macroscopic objects 

We now address the question raised in the beginning, i.e., on the validity of representing the 
quantum behavior of a macroscopic object by its center-of-mass dynamics, which we referred to as 
the 'COM axiom' for quantum dynamics of many body systems. To do this we consider a more 
general type of coupling between the system and the environment, e.g., coupling of the form f(xi)qj 
instead of XjQj, and examine if the COM variable dynamics separates from the reduced variable 
dynamics. 

For this purpose, let us note that if the function f(x) has the property J2iLi fi x i) = 
f(Xi) + g(X2, ...,Xn), for example f(x) = x or f(x) = x 2 , one can split the coupling between the 
system and environment into couplings containing the COM coordinate and the relative coordi- 
nates. Tracing out the environmental degrees of freedom qi, one can easily get the influence action 
which characterizes the effect of the environment on the system. 

However, the coarse-graining made by tracing out the environmental variables does not 
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necessarily lead to the separation of the COM and the relative coordinates in the effective action. 
When they are mixed up and can no longer be written as the sum of these two contributions, 
the form of the master equation will be radically altered as it would contain both the relative 
coordinate and the center-of-mass coordinate dynamics. One can work out how much of a change 
this would bring about in the COM dynamics, but at least we could say that when the coupling is 
not in these forms the COM axiom for macroscopic quantum dynamics no longer holds. 

Therefore we can conclude that for the N harmonic oscillators QBM model, the coupling between 
the system and the environment need be bi-linear, in the form Xtqj, for this axiom to hold. In that 
case, one can say that the quantum evolution of a macroscopic object in a general environment is 
completely described by the dynamics of the center-of-mass variable obeying a master equation of 
the HPZ type. 



C. Generalized standard quantum limit 

Another important issue of great interest to experimentalists is the generalized uncertainty 
relation of an N-body system at finite temperature. In our simple model the coupling between 
the system of NHO and the environment is only through the center-of-mass coordinate. If there is 
no mutual interaction between these N harmonic oscillators, only the center-of-mass coordinate is 
coupled to the environment and the remaining (N — 1) degrees of freedom are orthogonal. Hence 
the uncertainty function for the whole system will simply be Ujjz x (Ui) 1 ^' 1 , where Uhz is the 
uncertainty function for the 1HO QBM case given in [l^] and U± represents the quantum Heisenberg 
uncertainty relation for one pair of relative coordinate canonical variables Here we want 

to point out that in this situation, the effect of the environment, say the temperature of the heat 
bath, enters only through the center-of-mass coordinate in the form of Uhz, the remaining (iV — 1) 
pair of canonical conjugate variables are subject to quantum Heisenberg uncertainty relation U\. 
The time evolution and temperature dependence of the uncertainty function for an Ohmic bath 
was studied in great detail in [l?], Q] . 

If inter-particle interactions V%j 7^ exist among the N particles, then the uncertainty function 
Uhz x {Ui) N will be modified. Due to the interactions between the N harmonic oscillators, the 
uncertainty relation governing the (N — 1) pairs of relative coordinate variables (jQ,Pj) might 
be squeezed and rotated. However, if the number of harmonic oscillators N is large, and the 
interacting potential among the N-particles are short-ranged, or that the forces amongst them are 
very strong, so that the characteristic frequencies of these 'hard modes' are much higher than that 
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of the natural frequency of the COM modes, under these conditions, the quanta corresponding to 
the motion of the relative coordinates are not easily excited and the leading order contribution to 
the uncertainty function will be dominated by the center-of-mass degree of freedom. The details 
can be worked out from a perturbation analysis on the present results. (See, e.g., for treating 
one form of interaction.) 

D. Conclusion 

In this paper, we outlined the procedure to find a canonical transformation to transform from the 
individual coordinates (xj, Pi) to the collective coordinates (Xj, Pj), i = 1, N where X\,P\ are the 
center-of-mass coordinate and momentum respectively. We then proved that the potential Vij(xi — 
Xj) is independent of the center-of-mass coordinate X\. Then following the simple derivation 
of the master equation for 2HO in our previous work we showed that the system with variables 
(Xi,Pi) obeys a master equation of the same form as the HPZ equation. We gave the details 
of derivation of an exact non-Markovian master equation for the reduced density matrix of this 
system constructed with the heat bath variables integrated out. We stress that this result for the 
iV mutually interacting harmonic oscillators in a general environment is more than just a normal 
mode decomposition problem as in classical mechanics because there are interactions between the 
individual system oscillators and collective interaction with the quantum environment. 

This result is expected to be useful for the study of entanglement dynamics of multipartite 
particles and quantum to classical transition issues. Finally, we established a relation between the 
center-of-mass coordinate of N body harmonic oscillators and the well-known one oscillator QBM 
model. This provides a key step in establishing a microscopic theory for macroscopic quantum 
phenomena, a topic we intend to pursue further in the future. 

VI. APPENDIX 

A. Construction of the canonical transformation for general N . 

Given a system of ./V identical harmonic oscillators with equal mass M and intrinsic frequency 

n. 

i i 

Note that P { = M±i = and [xi,Pj] = iMij. 
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Goal: Find out {X;}, {Pi}, and {Mi} such that 

2M ^ 2 4 ^ 2M- ^2 * ' 

i i i 4 i 

where [Xj, Pj] = ihdij = [x i: Pj] and Xi = jj(xi + x 2 + ... + x/v), Pl = (Pi + P2 + ■■■ + Pn) and 
Mi = NM are the center-of-mass coordinate, total momentum, and total mass of the system 
respectively. Note that it's easy to see that [Xi,Pi] = ih. 



B. The explicit construction procedure: 

1. N = 2: 



Xi = i(xi+x 2 ), X 2 = xi-x 2 , 

Pi = Pi + P2, h = \(Pi- p 2 ), [x^ Pj] = ih£i 



(Mi,M 2 ) = (2M,y) 



Note that Pi = M { X { and Pi = Mii. 



fL + lL + l M tf X \ + l -Mtfx\ (40) 
2M 2M 2 2 1 y ' 

P " +^ + l(2M)n 2 Xl + ]-Atfxi (41) 



2(2M) 2(4f) 2 V ; 1 2 V 2 



2. X = 3 



: We can go from canonical variables (xi,x 2 ,X3) to (Xi = X cm , X2, X%) with the following 
transformation: Viewing {xi} and {Xj| as orthogonal base of the 3 dimensional vector space, our 
goal is to find a 3 x 3 linear transformation matrix T^, such that Xj = TijXj with det(Ty) = 1. A 
convenient choice is that since we have already construct the transformation of X = 2 case, we can 
just choose Xl = |(xi + x 2 + £3) and X 2 = x\ — x 2 . Since X3 must perpendicular with X\, X 2 , 
we have X3 = a 3 (xi + x 2 — 2x3). Since det(T) = 1, we then have 03 = \. 

Note that P3 = M3X3 hence proportional to X3, thus we have P3 = 63 (Pi + P 2 — 2P3). From 
[X3, P3] = ih, we have 63 = |. 

Xi = i(xi +x 2 + x 3 ) = X cm , X 2 = xi-x 2 , X 3 = i(xi +x 2 - 2x 3 ), 
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P 1 = (P 1 + P 2 + P 3 ), P 2 = i(Pi-P 2 ), P3 = ^(Pi + P 2 - 2P 3 ) 



and we have 



[Xi, Pj] = ihSij, [Xi, Xj] = [h, Pj] = 0. 

(Mi,M 2 ,M 3 ) = (2M,y,^ = ^M) 



A 2 



Pj Pf 1 
2M ' 2M + 2M + 2 

A 2 



Pi 

^ — h 



1 

2" 

l M, 



2(3 M) ■ 2(f) ' I + iW^ilyW^TW 



5. TV = 4; 



Note that in this case N=4=2 x 2. We can make use of our result of N=2. 
Since we require that X\ = \{x\ + x 2 + X3 + X4) in this case. We may first define 
yi = \{xi + x 2 ), y 2 = t;(x3 + X4), and 1/3 = (x\ - x 2 ),y 4 = (23 - x±). That means 
that we treat the whole 4 oscillators as two pairs of oscillators using the center-of- 
mass and relative coordinates of each pair. The corresponding conjugate momentum are 
Pyl = {Pi + P2), Py2 = (P3 + Pa),P V 3 = \{P\ ~ P2),P y 4 = \{Pz ~ Pa)- With this transformation, 
we have [yi,P y j] = ihdij. The effective masses for would be (2M, 2M, 4r, 4^). It's easy to see 
that 1/3 is perpendicular to all the other yi, and so is 2/4. Hence 1/3,1/4 are perpendicular to any 
linear combination of y\ and y 2 . That means by changing from Xi to yi we can decompose the 
original 4 dimensional vector space generated by Xi into direct sums of a two dimensional vector 
space generated by y\ and y 2 and two one dimensional vector spaces generated by y 3 and 1/4 
respectively. 

We can then make use of N = 2 result for y\,y 2 . Hence X\ = \{y\ + y 2 ) = 
\{x\ + x 2 + X3 + x 4 ),X 2 = yi — y 2 = \{x\ + x 2 - 23 - 24). The corresponding 

Pi = Pyl + Py2 = Pl+ P 2 + P3 + P±,P2 = \{Pyl ~ Py2) = \{Pl + ^2 " ^3 " Pi)- The 

effective masses are (2(2M), \{2M)) = {AM, M). 

Combing with the other two one dimensional vector spaces generated by y 3 and 2/4, we then 
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have the following transformation 



Xi = -{yi + y 2 ) = j(xi + x 2 + x 3 + x 4 ), 

X 2 = 2/i - 2/2 = ^{xi + x 2 -x 3 - x 4 ) 

^3 = V3=Xl- X 2 , 

X 4 = y 4 = x 3 - x 4 (42) 



with corresponding conjugate momentum 



P 1 = P yl + P y2 = P 1+ P 2 + P 3 + P 4 , 

^2 = - P y2 ) = \(Pi+P 2 -P 3 - Pa), 

Ps = P y3 = \{P 1 -P 2 ), 

P 4 = P yi = \{Pz-Pa) (43) 



and we have 



[Xi, P 3 ] = ih5 tJ , [X it Xj] = [Pi, Pj\ = 0. 

(Mi, M 2 , M 3 , M 4 ) = (AM, M, ™, ™ ) 

This procedure can be easily applied to any other cases where the number of oscillators are even. 
For example if N = 2xk. We can first define y\ = \(x\+x 2 ), y 2 = \(x 3 +x^), ...,y k = \{x 2 k-i+x 2 k) 
and y k+ i = (x 1 - x 2 ),y k+2 = (x 3 - x 4 ),...,y 2k = (x 2fc _i - x 2k ). The conjugate mo- 
mentum are P yl = (Pi + P 2 ),P y2 = (P 3 + P 4 ),...,P J/fe = (P 2k -i + P 2k ),P yk+ i = 
\(Pi - P 2 ),P yk+2 = \(P 3 - PA),--,P y2 k = \(P 2 k-\ ~ P 2 k) and the corresponding effective 
masses are (2M, 2M, ...2M, f, f, f). 

We then make use of the transformations for the N = k case to {yi}, i = 1, ...,k. By doing this 
we can have {Xi}, {Pi}, {M}, i = 1, k. Together with the Xj = yj, Pj = Pj,j = k + 1, 2k , 
we then have the complete transformation for N = 2k case. 

I N = 5: 

Because 5 = 4 + 1, we can make use of the result for N = 4 case. But note that in our 
construction we require that X\ = \(x\ + x 2 + x 3 + x 4 + X5) in this case, hence the corresponding 
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Pi = (P\ + P2 + P3 + P4 + P5) and effective mass is 5M. The corresponding X2,X3,Xt and 
P2, P3, P4 and effective masses are the same as those in the case of the N = 4 case. The remaining 
X5 has to be perpendicular to X^i = 1,2,3,4 and it's thus obvious that it must be a linear 
combination of the X[ = \(x\ + X2 + £3 + X4) in N = 4 case and X5 such that it is perpendicular to 
Xi = \(x\ + X2 + x% + X4 + X5). It's easy to see that X5 = a^(xi + X2 + X3 + X4 — 4xs). Hence the 
corresponding conjugate momentum is P5 = 65 (Pi + P2 + P3 + P4 — 4Ps). With the requirement 
that Xi = TijXj and det(Tjj) = 1, we have a$ = |. From the requirement [X5,Ps] = ih, we then 
have 65 = |. The corresponding effective mass is then given by M5 = = |M. 

Hence the transformation for N = 5 is given by 

Xi = i(xi + x 2 + x 3 + x 4 + x 5 ), 
X 2 = ^(xi + x 2 - x 3 - 24) 
X 3 = xi - x 2 , 
Xt = x 3 - x 4 , 

X 5 = ^(xi + X2 + x 3 + X4-4x 5 ), (44) 
with corresponding conjugate momentum 

Pi = Pl+P 2 +P 3 + P4 + P 5 , 
P 2 = ^(Pl+P 2 -P3-P 4 ), 

h = \{Pi-Pi), 

Pa = \(P*-Pa), 

P 5 = ^(Pi + P 2 + P 3 + P4-4P 5 ), (45) 

and we have 



[x^ Pj] = ihSij, [x^ Xj] = [h, Pj] = 0. 

M M AM 

(Mi,M 2 , M 3 , M 4 , M 5 ) = (5M, M, — , y , — ). 

This procedure can be easily applied to any other cases where the number of oscillators are odd. 
For example if N = 2k+l. We can make use of the results for N = 2k. Since we already have {xj} to 
{Xj}, i = I,..., 2k and these Xi are already orthogonal to each other. With the addition of X2fc+i we 
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only need to change X\ = 2 k+\ ( x i + x 2 + ■■■+x 2 k+x 2 k+i) and it's obvious that the X2k+i must have 
the structure X 2 k+i = a2fc+i(£i + x 2 + ••• + x 2 k — 2fc x 2 k+i)- With the requirement that Xi = TijXj 
and det(Ty) = 1, we have a 2k+1 = ^- Similarly P 2 fc+i = b 2k +i(Pi + P2 + ••• + Pik ~ 2kP 2k+1 ). 
The requirement [X 2k+ i, -P 2 fc+i] = ifo gives &2fc+i = 2T+1 nence the corresponding effective mass is 
M = hh±L M = = VM. 

We have thus provided a procedure to explicitly construct the canonical transformations for the 
canonical coordinates for any finite number of harmonic oscillators. The transformation for the 
effective masses are also given. With this canonical transformation, the structure of the Poisson 
brackets are preserved. 



C. Proof of the Lemma 



Lemma: If the potentials Vij among the system oscillators Xj and Xj are functions of X{ — Xj 
only, then = 0. Hence the total potential Yl^j-i^j Viji x i ~ x j) 1S independent of X\. 

Proof of the Lemma: The reason is very simple: regarding the original {xi},i = 1,...,N 
as an orthonormal coordinate basis which spans a iV-dimensional vector space. My construction 
for {Xi},i = 1,...,N is another set of complete orthogonal basis which contains the center-of- 
mass coordinate X cm = jj(x\ + ... + xn) as X\. In the {x{\ basis, X\ = -^(1, 1, 1) and 
Xi — Xj = (0, 0, 1, 0, .., —1, 0, ..,0) with the i-th entry equals to 1 and the j-th entry equals to 
-1. Hence it is obvious that Xi — Xj is orthogonal to X cm = X\. Hence x\ — Xj belongs to the 
(N — l)-dimensional vector space spanned by {Xj},j = 2, N of the new coordinates. Because 
{Xj}, i = 1, N is an orthogonal basis of the N-dimensional vector space V N , {Xj}, j = 2, N is 
an orthogonal basis for the (N-l)-dimensional subspace W N ~ l . Let L be the 1-dimensional vector 
space spanned by X x , then V N = L W N ~ 1 ,Ln W N ~ 1 = 0. Hence (a* - Xj) G W N ~ l and 
can be expressed in terms of {Xj},j = 2, N uniquely. Since we assume that the potentials are 
functions of only Xi — Xj, hence the potentials are functions of {X,}, j = 2, N only and are thus 
independent of X\. Q.E.D. 
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